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A recent paper (M. Rashdan, Phys. Rev. C 63, 044303 (2001)) introduces the new parameterization 
NL-RA1 of the relativistic mean-field model together with a new parameterization of the constant 
gap pairing model. Some conclusions of the paper may be doubtful as the pairing model is unrealistic 
and the known ground-state deformation of particular nuclei is neglected. 
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A recent study [1] introduces the parameterization NL- 
RA1 of the relativistic mean-field model (RMF). It is left 
open which data NL-RAl is fitted to, which prevents to 
relate its properties to the fitting strategy. NL-RAl is 
compared with early parameterizations as NL1 or NL- 
SH, more recent ones as NL3 or NL-Z2 are not consid- 
ered. Extrapolation of NL-RAl to superheavy nuclei con- 
tradicts earlier studies. As will be shown here, conclu- 
sions drawn in [1] may be doubtful as the pairing model 
used is unrealistic, and nuclei known to be deformed are 
calculated assuming spherical shapes. 

Ref. [1] employs two parameterizations of the constant- 
gap model, one of which (AI) has been used in many 
early applications of the RMF. The pairing matrix ele- 
ments are independent on the single-particle levels which 
is unrealistic for loosely-bound systems [2] as those dis- 
cussed in [1]. The pairing gap, related to the odd-even 
mass staggering, has to be parameterized as a function of 
N and Z [3]. Model AI describes the average behaviour 
of the pairing gap. Introducing about 30 parameters N c i 
and N C 2, model All attempts to incorporate the reduc- 
tion of the pairing gap around known magic numbers, but 
remains arbitrary for exotic systems where shells might 
be quenched or new shells occur. Model All misses the 
overall reduction of the pairing gap with A, see Fig. 1. 
For transactinides, gaps are overestimated by a factor of 
3. (There is currently much discussion about blocking 
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FIG. 1. Neutron pairing gaps A n from models AI and All 
in comparison with with experimental five-point gaps A' 5 '. 
The peak of A' 5 ' at closed shells is spurious, see [4]. The 
pattern of the gaps from model All around N ~ 180 is caused 
by the fact that the parameters AT c i and N C 2 of model All do 
not correspond to the nearest neutron shell closures. 
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FIG. 2. Ground-state quadrupole deformation 
I3 2 = <±n(r 2 Y 2 o)/CAArl) with r = 1.2 A 1/3 fm and relative 
error on the binding energy 5E of even-even transactinide 
nuclei. See [6] for a discussion of the uncertainty of SE. 

and mean-field contributions to calculated pairing gaps 
that are neglected here, see [4,5] and references therein. 
Those corrections cannot be easily incorporated into the 
simplistic constant gap model. Their contribution is usu- 
ally smaller than 20 % and decreases with A.) There is 
no justification for model All as it fails by construction 
to describe the size of the pairing gap, which is the key 
observable for pairing correlations. As the pairing gap 
determines the occupation of the single-particle states 
around the Fermi surface, most results presented in [1] 
are affected in one way or the other. 

In [1], the predictive power of NL-RAl for superheavy 
nuclei is tested for the heaviest known even-even nuclei as 
earlier done in [6] . Results in [1] differ significantly from 
[6]. One reason arc different pairing models, but there 
is a second one. In [1] it is not mentioned which shape 
degrees of freedom are accounted for. Repeating the cal- 
culations indicates that all results in [1] are obtained as- 
suming spherical shapes. This is consistent with exper- 
iment for Sn and Pb isotopes, but there is agreement 
among all successful mean-field models that the known 
superheavy nuclei are deformed. This is confirmed by 
experiment for selected isotopes up to 254 No [7]. In [1] 
the missing deformation energy (which is on the order of 
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FIG. 3. Single-particle spectrum of the protons, pairing 

gap A p in model All, difference A_E between the binding 
energy obtained with model AI and All, and two-proton gap 
<52 P (all in MeV) for N = 184 isotones. The uppermost panel 
also displays the occupation probability v\ for the protons 
in 114 calculated with model AIL The arrows denote the 
pairing gap A p , e p the Fermi energy. Results for 298 114 de- 
pend sensitively on the value of A p . Filled (open) markers de- 
note calculations where A p for 298 114 is calculated using the 
prescription for nuclei with N smaller (larger) than Z — 114. 
(Ref. [1] leaves it open which one to use. Results presented 
there correspond to filled markers). 

10 MeV or 0.5%) is replaced by the artificially increased 
pairing correlation energy from model AIL 

Using a more realistic state-dependent delta pairing 
force with parameters adjusted along the strategy of [4] 
and allowing for deformation change significantly the sys- 
tematics of 5E for transactinide nuclei, see Fig. 2. Com- 
paring with Fig. 17 in [1], all forces perform better. Sim- 
ilar changes can be expected for the values given in Fig. 
18 of [1] (see [8] for complications when calculating odd- 
A nuclei which are neglected in [1].) The change in SE 
when comparing NL1 with NL-Z and NL-Z2 reflects an 
improved center-of-mass correction [9] and the inclusion 
of data on exotic nuclei into an otherwise identical fit. 
NL-RA1 and NL-Z2, however, have the same good qual- 
ity for binding energies of transactinide nuclei in spite 
their very different nuclear matter properties. 

In the framework of mean-field models, a magic num- 
ber is associated with a large gap in the single-particle 
spectrum which causes a discontinuity in the system- 
atics of binding energies. Those (and other) disconti- 
nuities are filtered from data with the two-proton gap 
5 2p {N, Z) = E(N, Z - 2) - 2E(N, Z) + E(N, Z + 2) and 
the similar quantity for neutrons. 8 2p has to be taken 
with care, as it assumes the structure of the considered 
nuclei does not change, which is not necessarily fulfilled 
for heavy nuclei [10]. Ground-state deformation of some 



of the nuclei might quench 8 2p and has to be considered 
as soon as one wants to predict future data on these nu- 
clei. To demonstrate the non-existence of the spherical 
Z = 114 shell, however, spherical calculations are suffi- 
cient and enforce the validity of <52 P as a signature for 
magicity. Figure 3 displays the key quantities that reveal 
the origin of the large 8 2p for 298 114 found in [1]. Let 
us look at filled markers first. Single-particle energies 
in 298 114 do not change significantly when varying the 
pairing strength. The pairing gap for 298 114 from model 
All (A p = 1.7 MeV) is of the same order as the Z = 114 
gap in the single-particle spectrum (1.4 MeV). About 3 
protons occupy levels above the Z = 114 gap, inconsis- 
tent with the assumption of a major shell closure. By 
construction, A p drops by a factor 2 at Z = 114. This 
causes a discontinuity in the pairing correlation energy 
which is clearly visible when comparing binding energies 
obtained with models AI and AIL 

The discontinuity in the pairing correlation energy 
built into model All causes the large value for 52 P ( 298 114) 
found in [1], not the underlying shell structure. This is 
confirmed when moving the discontinuity of A p to the 
non-magic proton number Z = 112, see the open mark- 
ers in Fig. 3. S 2p is now peaked at 296 112, which has no 
closed spherical proton shell. <52 P cannot be used as a 
signature for shell closures when as the pairing gap and 
its fluctuations are of similar size as the spacing of single- 
particle energies. However, the size of the pairing gap in 
model All is unrealistic anyway. Calculations with more 
realistic pairing models do not show any significance for 
a major shell closure at Z = 114, consistent with [11]. 

The Z = 120 shell is not considered in the choice of 
N c i and N c2 in pairing model All, the pairing gap has 
huge mid-shell values there. This smears out the Z = 120 
shell effect found in [11] in terms of <52 P (c.f. Fig. 21 in 
[1]) but again does not affect the single-particle spectra. 
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